Absence of Asymptotic Freedom in Non-Abelian Models 



Adrian Patrascioiu 
Phystcs Department, University of Arizona, Tucson, AZ 85721, U.S.A. 



o 
o 
o 

(N 
D 
OO 



> ■ 

in 

O 

o : 
o . 

I 



X 



Erhard Seller 

Max- Planck- Institut fiir Physik (Wemer-Heisenberg-Institut) Fdhringer Ring 6, 80805 Munich, Germany 

(February 1, 2008) 

The percolation properties of equatorial strips of the two dimensional 0(3) nonlinear o model 
are investigated numerically. Convincing evidence is found that a sufBcently narrow strip does not 
percolate at arbitrarily low temperatures. Rigorous arguments are used to show that this result 
implies both the presence of a massless phase at low temperature and lack of asymptotic freedom 
in the massive continuum limit. A heuristic estimate of the transition temperature is given which 
is consistent with the numerical data. 
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One of the crucial unanswered problems in particle and 
condensed matter physics concerns the phase diagram of 
the two dimensional {2D) nonlinear a models m. It is 
widely believed that the nonabellan models with iV > 3 
are in a massive phase for any finite inverse temperature 
/3 and that this is intimately related to their perturba- 
tive asymptotic freedom. Over the years we have brought 
forth many reasons why we think that these beliefs are 
unfounded |^ ^ , but the absence of a mathematical proof 
combined with ambiguous numerical results left the is- 
sues wide open. In the present letter we would like to 
provide convincing numerical evidence that in fact the 
2D 0(3) model possesses a massless phase for sufficiently 
large (3 and give a rigorous proof that this is incompat- 
ible with asymptotic freedom in the massive phase. We 
will also give a heuristic explanation of why and where 
the phase transition happens. 

The models we are considering consist of unit length 
spins s taking values on the sphere S^^^, placed at the 
sites of a 2D regular lattice. These spins interact ferro- 
magnetically with their nearest neighbors. Let {ij) de- 
note a pair of neighboring sites. We will consider two 
types of interactions between neighbouring spins: 

• Standard action (s.a.): Hij = —s{i) ■ s{j) 



• Constrained action (c.a.): Hij = ~s{i) ■ s{j) for 
s{i) • s{j) > c and Hij ^ oo for s(i) • s{j) < c for 
some c £ [—1, 1). 

Almost a decade ago we showed ||] that one can 
rephrase the issue regarding the existence of a soft phase 
in these models as a percolation problem and in fact this 
is the reason we introduced the c.a. model (please note 
that the c.a. model has the same perturbative expansion 
as the s.a. model and possesses instantons). Namely let 
e — \/2{l — c) and Se the set of sites such that \s-n\ < e/2 
for some given unit vector n. Our rigorous result was that 
if on the triangular (T) lattice the set Se does not con- 
tain a percolating cluster, then the 0{N) model must be 



massless at that c. For the abelian 0{2) model we could 
prove the absence of percolation of this equatorial set Se 
for c sufficiently large (modulo certain technical as- 
sumptions which were later eliminated by M. Aizenman 
Q). For the nonabellan cases we could not give a rig- 
orous proof. We did however present certain arguments 
[^,|| explaining why the percolating scenario seemed un- 
likely. 

In this letter we will give convincing numerical evi- 
dence that a sufficiently narrow equatorial strip does not 
percolate for any c. We will also show that via a rigorous 
inequality derived by us in the past the existence of 
a finite (3crt in the s.a. model on the square (S) lattice is 
incompatible with the presence of asymptotic freedom in 
the massive continuum limit. 

For clarity we will review first a few crucial points re- 
garding percolation in 2D: 

a) Let A be a subset of the lattice defined by the spin 
lying in some subset A C S^~^ and A its complement. 
Then with probability 1 A and A do not percolate at the 
same time. (For this point it is crucial that the lattice 
is self matching, hence our use of the T instead of the S 
lattice; on the latter an ordinarily connected cluster can 
be stopped by a cluster connected only Tkr-wise, i.e. also 
diagonally.) This fact has been proven rigorously only 
for special cases like the -|- and — clusters of the Ising 
model, but is believed to hold generally. 

b) If neither A nor its complement A percolate, then the 
expected size of the cluster of A attached to the origin, 
denoted by (A) of A, diverges; the same holds for its 
complement A (Russo's lemma [|l^). (If A percolates, 
then (A) is finite.) 

Let then Ve (the union of the polar caps) be the com- 
plement of Se (the equatorial strip of width e) . According 
to the discussion above, either Se percolates, or Pe perco- 
lates, or neither Se nor Pe percolates and then both have 
divergent mean size (we shall call this third possibility in 
short formation of rings) . Consider a sequence of tori of 
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increasing size L and the mean cluster size of the set A 
corresponding to a subset A C 5"^"^ of positive measure. 
If A percolates {A) — 0{L^), if its complement percolates 
(A) will approach a finite nonzero value, and if A forms 
rings (A) ~ 0{L^~^) for some t] > 0. Therefore, if we 
define the ratio 



(1) 



for L — > cxD it should either go to 0, or to oo or to some 
finite nonzero value depending on which one of the three 
possibilities is realized; the latter possibility assumes that 
T] is the same for and Se, as indicated by our numerics 
and consistent with scaling. 

In Figjl] we show the numerical value of the ratio r as 
function of c for /3 = for four values of e for the c.a. 
model on a T lattice. 
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FIG. 1. Ratio {P^)/{Se) for various e values versus c 



The results were obtained 
investigation using an 0(3) 
Wang cluster algorithm and 
20,000 lattice configurations 
ments. For each value of e 
and 640 (for e = .78 we also 
distinct regimes are manifest 
of e investigated: 



from a Monte Carlo (MC) 
version of the Swendsen- 
consist of a minimum of 
used for taking measure- 
we studied L = 160, 320 
studied L = 1280). Three 
for each of the four values 



• For small c, r is increasing with L, presumably di- 
verging to oo (region 1). 

• For intermediate c, r is decreasing with L, presum- 
ably converging to (region 2). 

• For c sufficiently large depending upon e, r becomes 
independent of L, just as it does at the crossings 
from region 1 into 2. 

Consequently for these values of e for small c perco- 
lates, for intermediate c S^, percolates and for sufficiently 
large c both P^ and S'e form rings. In Fig.^ we present 
the phase diagram of the c.a. model on the T lattice for 
P = 0. The dashed line D represents the minimal equa- 
torial width above which percolates. For c = — 1 (no 



constraint) the model reduces to independent site per- 
colation, for which the percolation threshold is known 
rigorously to be e = 1. The rest of the diagram repre- 
sents qualitatively the results of our investigation of the 
ratio r, such as those shown in Fig.^ Two features of 
this diagram are worth emphasizing: 

1. An equatorial strip of width less than approximately 
e = .76 never percolates. 

2. For approximately c > 0.4 a new phase opens up in 
which both P^ and Se form rings (the dotted line sepa- 
rates it from percolation of Pg). 




FIG. 2. Phase diagram of the 0(3) model on the T lattice 

For clarity let us briefly review the argument indicat- 
ing that this phase diagram is incompatible with the exis- 
tence of a massive phase for arbitrarily large c. Choosing 
an arbitrary unit vector n, one introduces Ising variables 
a = ±1. The s.a. Hamiltonian becomes: 
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where S||(i) = s{i) ■ n and sj_{i) = (s(*) x n) x n. 
One thus obtains an induced Ising model for which the 
Fortuin-Kastleyn (FK) representation |^ is applicable. 
In this representation the Ising system is mapped into 
a bond percolation problem: between any like neigh- 
boring Ising spins a bond is placed with probability 
p — 1 — exp(— 2/3s||(z)s||(j)). For the c.a. model a bond 
is also placed if after flipping one of the two neighboring 
Ising spins the constraint s{i) ■ s{j) > c is violated. The 
FK representation relates the mean cluster size of the 
site clusters joined by occupied bonds (FK-clusters) to 
the Ising magnetic susceptibility. In a massive phase the 
latter must remain flnite. Hence, if the FK-clusters have 
divergent mean size, the original 0(3) ferromagnet must 
be massless (the Ising variables a are local functions of 
the originally spin variables s). 

Now notice that by construction for the c.a. model 
the FK-clusters with say a = +1 must contain all sites 
with s{i) -n > -^(1 — c)/2. Therefore the model must be 
massless if clusters obeying this condition have divergent 
mean size. But the polar set P^ consists of two disjoint 
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components P+ (north) and P~ (south). For c > (1 — 
e^)/2 there are no clusters containing both elements of 
P+ and P~ . Hence if for such values of c clusters of 
form rings, so do clusters of P+ separately and hence 
the 0(3) model must be massless. The curve C given by 
c = (1 - e2)/2 is the solid hnc in FigJ|. The point T at 
the intersection of the curves D and C gives an upper 
bound for Ccrt, the value of c above which the c.a. model 
must be massless. 
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FIG. 3. The ratio of the mean cluster size of a polar cap of 
height .75 to that of an equatorial strip of the same height 

To verify the phase diagram in Fig.|^ we also measured 
(at f3 ~ 0) the ratio of the mean cluster size of the set 
P,t with e' = .5 to that of the set Se with e = 0.75, for 
which the equatorial strip never percolates (Fig.^ (e and 
e' ar chosen such that the two sets have equal density). 
The data indicate that for some intermediate values of c 
P^ forms rings while S'^ has finite mean size; this region 
terminates around c = 0.4, where also starts forming 
rings. The larger average cluster size of the polar cap 
compared to the strip of the same area is probably due 
to the fact that the strip has a larger boundary than the 
polar cap. This is in agreement with a general conjecture 
stated in jTj, namely that for c sufficiently large, if two 
sets have equal area but different perimeters, the one with 
the smaller perimeter will have larger average cluster size. 
For the case at hand, this is apparently true for all values 
of c. 

The general belief, which we criticized in ref. §,is 
that there is a fundamental difference between abelian 
and nonabelian models. To test this belief we studied 
the ratio r for the c.a. 0{2) model on the T lattice. 
The phase diagram is shown in Fig.^. Since in the 0{2) 
model the set P e can also be regarded as a set Si where 
e = ■\/(4 — e^), certain features of that diagram follow 
from rigorous arguments. For instance it is clear that 
in the c.a. model there exist two curves C and C and 
in the region to their right the model must be massless 
The precise location of the curves D (or D) must 
be determined numerically, something which we did not 
do. We did verify though that the ring formation region 



begins around c = —0.5. 

In our opinion the arguments and numerical evidence 
provided so far give strong indications that the c.a. mod- 
els on the T lattice possess a massless phase. Universality 
would suggest that a similar situation must exist for the 
s.a. models on the T and S lattices. To test universality 
we measured on the S lattice the renormalized coupling 
both on thermodynamic lattices in the massive phase and 
in finite volume in the presumed critical regime (as in 
Ip^). Our data for the c.a. model on the S lattice only 
determine an interval (about .5 to .7) in which the mass- 
less phase of the model sets in; we tried to see if we could 
get a similar L dependence for the renormalized coupling 
in the s.a. model at a suitable /? as for c = .61 in the c.a. 
model at /3 = 0. This seems to be indeed the case for 
f3 roughly 3.4. We went only up to L = 640, hence this 
equivalence between c and P should be considered only as 
a rough approximation, but there seems to be no doubt 
that the two models have the same continuum limit. 




FIG. 4. The phase diagram for the 0(2) model on the T 
lattice 

It is intersting to note that there is a heuristic expla- 
nation for both the existence of a massless phase in the 
s.a. 0(3) model and for the value of /3crt- Indeed it 
is known rigorously that in 2D a continuous symmetry 
cannot be broken at any finite f3. In a previuos paper y] 
we showed that the dominant configurations at large /3 
are not instantons but superinstantons (s.i.). In principle 
both instantons and s.i. could enforce the 0(3) symme- 
try. In a box of diameter R the former have a minimal 
energy Emst = 47r |jl^ while the latter Esj. = S^n/lnR, 
where 5 is the angle by which the spin has rotated over 
the distance R. Now suppose that (3crt is sufficiently 
large for classical configurations to be dominant. Then 
let us choose 6 — 27r (restoration of symmetry) and ask 
how large must R be so that the superinstanton config- 
uration has the same energy as one instanton. One finds 
In R = TT^. But in the Gaussian approximation 

(s(0) • s(a;)) « 1 - -!-lnx (3) 
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Thus restoration of symmetry occurs for Inx tt/?. 
This simpleminded argument suggests that for (3 > tt 
instantons become less important than s.i.. Now in a 
gas of s.i. the image of any small patch of the sphere 
forms rings, hence the system is massless. While this 
is not a quantitative argument, we believe it captures 
qualitatively what happens: a transition from localized 
defects (instantons) to s.i.. 

Next let us discuss the connection between a finite Pert 
and the absence of asymptotic freedom. It follows from 
our earlier work concerning the conformal properties of 
the critical 0(2) model We refer the reader for de- 
tails to that paper and give only an outline of the argu- 
ment. The s.a. lattice 0{N) model possesses a conserved 
isospin current. This currrent can be decomposed into a 
transverse and longitudinal part. Let F'^{p) and F'"{p) 
denote the thermodynamic values of the 2-point functions 
of the transverse and longitudinal parts at momentum p, 
respectively. Using reflection positivity and a Ward iden- 
tity we proved that in the massive continuum limit the 
following inequalities must hold for any p ^ 0: 

< F'^ip) < F'^iO) = F^{0) < F^{p) = 2I3EIN (4) 

Here E is the expectation value of the energy 

E = {s{i) ■ 

at inverse temperature /3. Since i? < 1 it foUws that if 
Pert < oo F'^{0) — F'^{p) cannot diverge for p — > oo as 
required by perturbative asymptotic freedom. In fact, 
for Pert = 3.4 (which is a plausible guess) F'^{p) must be 
less than 2.27, in violation of the form factor computation 
giving f ^(0) - F^{oo) > 3.651 ll4|. 

Since the implications of our result, that for the c.a. 
model a suffciently narrow equatorial strip never perco- 
lates, are so dramatic, the reader may wonder how cred- 
ible are the numerics. The usual suspect, the random 
number generator, should not be important since preci- 
sion is not the issue here. The only debatable point is 
whether our results represent the true thermodynamic 
behaviour for L — > oo or are merely small volume arte- 
facts. While we cannot rule out rigorously that scenario, 
certain features of the data make it highly implausible: 

• Small volume effects should set in gradually, while 
the data in Fig.^ indicate a rather abrupt change 
from a region where r is decreasing with L to one 
where r is essentially independent of L. 

• For c ^ 1 at fixed L, r must approach the 'geo- 
metric' value r = 2/e — 1. As can be seen, in all 
the cases studied, throughout the 'ring' region r is 
clearly larger than this value, while it should go to 
if percolated. 

• In Fig.| there is no indication of the ratio going to 
for increasing L. Moreover the dramatic change 
in slope around c = .4 indicates that the polar cap 
Pe' starts forming rings at a smaller value of c than 
the equatorial strip Se- 



Thus we doubt very much that the effects we are see- 
ing represent small volume artefacts. Moreover, if Sz 
remained massive at low temperature and in fact an ar- 
bitrarily narrow equatorial strip percolated, one would 
have to explain away our old paradox : if such a nar- 
row strip percolated, an even larger strip would percolate 
and on it one would have an induced 0(2) model in its 
massless phase, in contradiction to the Mermin- Wagner 
theorem. 

Consequently we are confident that the phase diagram 
in Fig.|^ represents the truth, that a soft phase exists 
both in the s.a. and the c.a. model and that the massive 
continuum limit of the 0(3) model is not asymptotically 
free. In a previous paper 0] we have already shown that 
in nonabelian models even at short distances perturba- 
tion theory produces ambiguous answers. The present 
result sharpens that result by eliminating the possibility 
of asymptotic freedom. 
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